Daisy models: Semi-Poisson statistics and beyond 
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Semi-Poisson statistics are shown to be obtained by remov- 
ing every other number from a random sequence. Retaining 
every (r + l)th level we obtain a family of sequences which we 
call daisy models. Their statistical properties coincide with 
those of Bogomolny's nearest-neighbour interaction Coulomb 
gas if the inverse temperature coincides with the integer r. In 
particular the case r = 2 reproduces closely the statistics of 
quasi-optimal solutions of the traveling salesman problem. 

PACS numbers: 03.65. Sq, 05.45.+b 



The transition from order to chaos in a classical system 
is generically reflected by a transition in spectral statis- 
tics for the corresponding quantum system from Poisson 
to GOE statistics. Such a transition was first considered 
by Porter and Rosenzweig Q in a model that features a 
parameter that depresses the off-diagonal elements of a 
GOE until they are zero and we have Poisson statistics. 
This model is amenable to analytical treatment but 
does not reflect the properties of dynamical systems very 
well. Band matrices were introduced later || and have 
been quite successful in the description of many situa- 
tions ||. A semi-classical ansatz by Berry and Robnik 
H was shown to work very well if applied to the long- 
range behaviour of the two-point function || . 

More recently a different kind of transitional behaviour 
was discovered near derealization transitions [Q and in 
pseudo-integrable systems P,pl; ^ ^ s commonly known 
as semi-Poisson statistics. Bogomolny, Gerland and 
Schmit developed a level dynamics for this type of spec- 
tra by limiting the usual 1 — D Coulomb gas model to 
nearest-neighbour interactions and considering an inverse 
temperature (3 — 1 ||. They also pointed out that 
the nearest-neighbour spacing distribution can be repro- 
duced by an interpolation procedure in a Poisson spec- 
trum 

The first purpose of the present note is to show that the 
statistics of a semi-Poisson spectrum are reproduced ex- 
actly by a model where every other level is dropped from 
a Poisson spectrum. The idea for such a model derives 
from well-known results that relate the superposition of 
two GOEs 
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to a GUE or one GOE to a GSE by the same proce- 
dure. We shall term such models daisy models. As in the 
above cases a dynamical link is not established. Rather 
we calculate properties of the spectra and find that they 
coincide. This is certainly no dynamical explanation of 
the properties of the physical systems that display semi- 
Poisson statistics, but neither is any of the level-dynamics 
models. The main advantage is that it is a simple model 
for which it is very easy to compute any statistic. 

It is further interesting to note that this procedure does 
not only yield one new type of statistical spectra but 
rather an entire family of daisy models of rank r where 
r is the number of levels dropped between each retained 
level. Their statistical properties can be easily calculated 
for any r. As we shall see, the dependence on rank is 
exactly the same as on the inverse temperature (3 in the 
nearest- neighbour interaction Coulomb gas. Thus we find 
that all integer values of this parameter correspond to a 
daisy model. 

No link of such statistical spectra to quantized dynam- 
ical systems is known for r > 1, but the case r = 2 dis- 
plays surprising similarity to the statistical distribution 
of distances between cities along a quasi-optimal path of 
the traveling salesman problem ]l0| . It may be worth- 
while to note the relation of this problem to spin glasses 
p"l[ though we shall not discuss this aspect. 

The semi-Poisson spectra display a nearest-neighbour 
spacing distribution 



P(s) = 4sexp(-2s) 



(1) 



and a long-range behaviour of the two-point function de- 
fined by a number variance 



E 2 (L) = L/2. 



(2) 
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The nearest-neighbour interaction 1 — D Coulomb gas 
model is defined as follows: We have A^ + 2 particles with 
positions Xj in an interval of size L with the interaction 

V(x , xi,X2, xn+i) = ~y] log(a; 4 - Xi-i) (3) 



and the condition = xq < x\ < .... < xn < Xjv+i = 
L. This model has the following nth- neighbour spacing 
distribution for inverse temperature (3 

(R + 1 \{P+ 1 ) n 
S) = T([P+l]n) s{P+1)n ^ ^PH/ 3 + !H ( 4 ) 
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In particular for (3 = 1 and n = 1 we obtain the nearest- 
neighbour distribution for the semi-Poisson (eq. 1). 

For the daisy model of order r as defined above we ob- 
tain the nth-neighbour spacing distribution by rescaling 
the rnth-neighbour distribution of a random sequence as 

(r -A- 1 )( r + 1 ) n 

P r (n,s) = [ J >— ^'"^expKr+Di (5) 
T([r + l\n) 

The same rescaling argument yields for the asymptotic 
behaviour with L ^> 1 of the number variance for rank r 
model, Sr(^) ~ L/(r + 1). Using the expression (4.41) 
of ref. Q for the two-point function we obtain for the 
number variance 




£?(£) = £• 



2L 



1 U (1 - Wj) ( r + 1) 



(6) FIG. 1. Fig. 1: The number variance E for the traveling 

salesman problem (dashed line) compared to the one of the 
rank 2 daisy model (full line). 
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■exp([Wj-l}{r+l)L) 



where Wj — exp(-^rf ) are the r + 1 roots of unity. All 
sums in this expression are real and the fist two can be 
summed as follows: 

Consider the function P(x) = p= = IE (x — 

Wj), with Wj as above, and its logarimic deriva- 
tives. We note that \og(P(x))] = r/2 = 

£,\og{P{x)) 



-r +ir 
12 



E^iTzkvJ and 

^y~j=i (i-lv ) 2 • Hence, it is possible to write for the first 
sum T r 3 Wj - V W '~ 1 + V r 1 - -r/2 
and for the second one 



\ " 



'Wj-1 



4rf (1-W^i) 2 f-f (1-W,-) 2 ^(1- 



r 1 
^ (1 - Wi 



(7) 



r 

= ~2 
We thus obtain 
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-r(r + 2) 
12 ' 



L r(r + 2) 



+ 1 6(r + 1)' 



(r + 



l) 2 tt(l-^ eXP 



([^■-l](r + l)i), 



which corroborates the asymptotic behaviour for L>1. 
In particular, for r = 2 we have 



L 4 
3 + 27 



1 r 3 ^n / 9 n 
1 -cos(— Z)exp(--L) 



which is depicted by the solid line in fig 1. 

If we compare equations (4) and (5) we find that they 
coincide for (3 — r. Thus, we link the inverse temperature 
in the Coulomb gas model to the number of discarded 
elements in a rank r daisy model starting from a random 
sequence. The fact that semi-Poisson statistics can be 
interpreted in this simple way is of interest because the 
properties of daisy models over random sequences are 
easy to calculate and because it may shed some light 
into the 1 — D Coulomb gas dynamics. 

On the other hand, we may ask if models of rank larger 
than 1 are relevant. Some of us pointed out recently 
|Tof that the statistical distribution of distances between 
cities along a quasi-optimal path of the traveling sales- 
man problem displays characteristic features that can be 
analyzed with the tools of spectral statistics, but cannot 
be understood in terms of any of the usual random ma- 
trix models. In particular, the long-range behaviour of 
the number variance E 2 and the correlation coefficient 
between adjacent spacings seem quite incompatible with 
band matrices or Porter- Rosenzweig type models [ [To] . 
We shall therefore investigate whether we obtain a bet- 
ter agreement with a daisy model. 

The data for the traveling salesman problem are ob- 
tained for an ensemble of 500 maps of 500 cities using 
simulated anneling. We first compared with the pseudo- 

(8) Poisson model and found no agreement. But the slope of 
the number variance suggests that we should rather com- 
pare with the rank 2 daisy model for a random sequence. 

In figure 1 we compare the number variance with eq. 
(9) for the rank 2 model and we find a remarkable agree- 
ment. If we fit the asymptotic of the number variance 
E 2 (^) with a straight line for the interval 2.0 < L < 9.0 
we obtain a slope of 0.358 ± 0.001 near to the value 1/3. 
Furthermore, the ith-neighbour spacing width <jj for the 
rank 2 model is equal to a rescaled 2zth-neighbour distri- 

(9) bution of a Poisson ensemble. In fig 2 we compare these 
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FIG. 2. Fig. 2: The width cr, of ith-neighbour spacing 
distribution of the traveling salesman problem (crosses) com- 
pared to that of the rank 2 daisy model (full line) . 



widths with the ones obtained for the traveling salesman 
problem in fic[| and find similar agreement. For the cor- 
relation coefficient in the rank 2 model we expect zero, 
which is quite near to the value .036 ± .002 of ref. |jTo| . 
We also compared the properties of skewness and curtosis 
used commonly to detect properties of the three and four 
point functions. The results for the rank 2 model were 
obtained numerically and the comparison is displayed in 
figs 3a and 3b. Similarly the nearest-neighbour spac- 
ing distributions are compared in fig. 4. The agreement 
is certainly comparable to the one obtained for pseudo- 
integrable systems with the semi-Poisson statistics. 

Summarizing, we have shown that semi-Poisson statis- 
tics can be obtained from a very simple model without 
any dynamical implications. This model actually per- 
tains to a family of models that seem to be relevant in 
situations where the usual banded matrix models and the 
Porter Rosenzweig model are grossly inadequate. We ob- 
tain this family by retaining every (r + l)th level of a 
random sequence. A similar selection process could be 
performed for the classical ensembles of Cartan (e.g. the 
GOE) p2| . Whether this leads to useful results beyond 
the two cases mentioned above is an open question. 
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FIG. 3. The skewness (a) and the curtosis (b) of the 
traveling salesman problem (crosses) compared to those of 
the rank 2 daisy model (diamonds). 
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